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Abstract. - This is an account of the computation of X-ray microstrain in a polycrystal with 
anisotropic elasticity under uniaxial external load. The results have been published in the article 
"Microstrain in nanocrystalline solids under load by virtual diffraction", at Europhysics Letters 
89, 66002 (2010). The present information was submitted to Europhysics Letters as part of 
the manuscript package, and was available to the reviewers who recommended the paper for 
publication. 



Introduction. — This text accompanies a recent pub- 
£h location [I] in which the present authors propose a method 



for computing virtual x-ray diffractograms for nanocys- 



i— talline materials under uniaxial load based on molecular 
dynamics simulation data. While the increase in diffrac- 



tion microstrain during deformation is generally taken as 



evidence for the generation of lattice dislocations, the vir- 
r^j- tual diffraction data in Ref. [1] show extra microstrain even 
0^ at small load, before the onset of lattice dislocation activ- 
\l ity. Reference [T] argues that the microstrain data have a 
CO natural explanation in the elastic response of a heteroge- 
^Dijeous medium to uniaxial external load. The results im- 
^^ ply that the conclusions of previous experimental in-situ 
. . diffraction data for nanocystalline metals may deserve a 
^ critical examination. 

, Here, we account for the computation of X-ray micros- 

train in a polycrystal with anisotropic elasticity under uni- 
5^ axial external load from the perspective of the mechanics 

of heterogeneous elastic media. 

General approach to microstrain. We consider a 
space-filling array of grains with random crystallographic 
orientation, in a way that the polycrystalline material 
is macroscopically isotropic, and we restrict attention to 
purely elastic deformation. We allow for anisotropic elas- 
tic response of the crystal lattice of the individual grains, 
consistent with cubic symmetry. In a diffraction experi- 
ment the heterogeneous strain in the polycrystal will lead 
to broadening of the Bragg reflections which is measured 



by the diffraction microstrain, as described below. 

Our analysis uses crystal coordinates. Consider diffrac- 
tion at a Bragg reflection with Miller indices hkl, and let 
the vectors u*hkh Vhki and Whki form an orthonormal ba- 
sis set with Uhki = (h, k, l)/^/h 2 + k 2 + l 2 the unit vec- 
tor along the zone axis that is aligned with the scat- 
tering vector, q. For instance, for the (422)-reflection 
one may take u = (2, 1, l)/v6, v = (0,— l,l)/v2 and 
w = (—1,1, 1)/a/3- Since the diffraction experiment has q 
normal to the load axis, a general expression for the unit 
vector n along the load axis is 



fihki (if) = v hk i cos <p + w hk i sin (p 



(1) 



with ip an orientation variable. In other words, at any 
given Bragg reflection hkl the crystals which satisfy the 
Bragg condition are responding to an external load that 
can take on different directions as seen within a crystal 
coordinate system. When ip varies between and 2tt, then 
the vector fihkiif) covers the entire domain of relevant 
orientations of the external load axis. The uniaxial stress 
due to the respective load is given by the tensor 



T hk i{p) = cr ri hk i(tp) ® n hk i(tp) 



(2) 



Here, a denotes the magnitude of the stress. 

The interference function of the polycrystalline aggre- 
gate is the superposition of Bragg reflections of crystals 
with a variety of load directions in crystal coordinates, 
as discussed above. The Bragg reflection position of any 
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given crystallite is determined by the lattice parameter, 
measured in the direction of q. Since the strain - in other 
words, the relative change of the lattice parameter as com- 
pared to a stress-free reference configuration - depends on 
the load direction, there will generally be a distribution of 
lattice parameters and, hence, a distribution of reflection 
positions. It is this fact that links the reflection broaden- 
ing to the load via the anisotropic elastic response of the 
crystal lattice. 

Specifically, the reflection broadening is measured by 
the microstrain, ehki, which depends on the projection, 
e/ife/(v)i of the strain E onto Uhki via 

I /" 27r 

e lki = ^J o (thki{<p) - (e hk i)) dip (3) 

In the following we compute ef lkl using first the Rcuss and 
then the Kroner approximation. 

Reuss approximation. — As a first step towards 
analysing the strain, consider the Reuss approximation, 
where the stress is assumed uniform throughout the poly- 
crystal. The local stress in any grain is then given by 
equation (J2J), and the strain is 

E% kl (<p)=S:T hkl (<p) (4) 

with S the compliance tensor. The superscript refers to 
the Rcuss aproximation. The strain projected on the scat- 
tering vector is here 

ehki((f) = u hkl ■ F,f kl ((p) ■ u hkl . (5) 

By using equations (H]) and along with tabulated 
values for the single-crystal compliance coefficients (rep- 
resented, for instance, by Cn, C12, C44 for a cubic crystal 
lattice), equations (|3|) for the square of the microstrain 
value ehki is readily evaluated at each Bragg reflection. 

Kroner approximation. To evaluate the micros- 
train in the Kroner approximation [2 , 3 we take the grain 
of orientation h,k,l,ip as an inclusion that is embed- 
ded in an elastically isotropic continuum with the self- 
consistently averaged (Kroner-) elastic constants, repre- 
sented here by S. The entries of the tensor S can be com- 
puted by means of equations (21) and (22) in reference [2J. 
Sufficiently far from the inclusion, the continuum is then 
strained by 

iWv) = S : T hH (cp). (6) 

The strain within the grain is uniform, and it is given by 
(cf. equation (14) in reference [5]) 

Eww(v)= (s+twwfeO) :T hkl (cp). (7) 

The quantity th k i((p) describes an excess compliance due 
to the interaction between grain and matrix. It can 
be computed from the single-crystal elastic constants by 
means of equations (18) and (19) in reference [2J. 



Table 1: The dimensionless microstrain coupling coefficients 
b for Pd (Sydow potential) in different crystallographic ori- 
entations (hkl). Values derived from Reuss (b R ) and Kroner 
(6 K ) models are distinguished by superscripts. The bottom line 
shows the root-mean square of the coefficients. 
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(110) 


0.213 


0.106 


(111) 








(210) 


0.137 


0.068 


(211) 


0.071 


0.035 


(221) 


0.126 


0.063 


(310) 


0.077 


0.038 


(311) 


0.056 


0.038 


(321) 


0.130 


0.065 


(331) 


0.170 


0.085 


RMS 


0.118 


0.059 



Analogously to equation (O, the strain projected on the 
scattering vector is here 

ehkliv) = u hk i ■ EwmOp) ' UhM • (8) 

By using equations along with the single-crystal 

compliance coefficients, we can again evaluate equation J3]) 
for the square of the microstrain value ehki at each Bragg 
reflection. 

Results. — The microstrain values ehki computed 
from our model arc linear in the stress magnitude, a. We 
can therefore specify a stress-independent microstrain pa- 
rameter, b, defined so that 

ehki = bhkia/Y . (9) 

Here Y is taken as the macroscopic Young modulus of the 
respective model, Reuss or Kroner. 

Table [JJ shows the values obtained with the Reuss and 
Kroner models. Values for the single crystal stiffness 
were those of the Sydow potential for Pd [J], namely 
C11 = 292.7 GPa, C 12 = 188.9 GPa, C 44 = 125.0 GPa. 
The corresponding Kroner elastic constants were obtained 
as Y = 235 GPa, G = 88.7 GPa, v = 0.325 where V, G, v 
denote the Young modulus, the shear modulus, and the 
Poisson ratio, respectively. . 
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